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1 Introduction 

High Energy collisions of nuclei usually lead to the emission of evaporation 
neutrons from both "beam" and "target" nuclei. At the RHIC heavy ion col- 
lider with lOOGeV/u beam energy, evaporation neutrons diverge by less than 

2 milliradians from the beam axis. Neutral beam fragments can be detected 
downstream of RHIC ion collisions (and a large aperture Accelerator dipole 
magnet) if ^ < 4 mr but charged fragments in the same angular range are 
usually too close to the beam trajectory. In this 'zero degree' region produced 
particles and other secondaries deposit negligible energy when compared with 
that of beam fragmentation neutrons. 

The purpose of the RHIC zero degree calorimeters (ZDC's) is to detect neu- 
trons emitted within this cone along both beam directions and measure their 
total energy (from which we calculate multiplicity). The ZDC coincidence of 
the 2 beam directions is a minimal bias selection of heavy ion collisions. This 
makes it useful as an event trigger and a luminosity monitor [1] and for this 
reason we built identical detectors for all 4 RHIC experiments. The neutron 
multiplicity is also known to be correlated with event geometry [2] and will 
be used to measure collision centrality in mutual beam interactions. 



2 Design Goals 

The RHIC ZDC's are hadron calorimeters. Their longitudinal segmentation 
(2 A/, 50Xo) is determined by practical, mechanical considerations. Electro- 
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Fig. 1. Plan view of the collision region and (section A- A) "beam's eye" view 
of the zdc location indicating deflection of protons and charged fragments ( with 
Z/A~ 1 downstream of the "DX" Dipole magnet. 

magnetic energy emission into this region is predicted to be negligible so this 
measurement is not emphasized in our design. Since the spatial distribution of 
neutrons emitted in the fragmentation region carries only limited information 
about the collision, the calorimeters are built without transverse segmentation. 

The Forward Energy resolution goal was determined by the need to cleanly 
resolve the single neutron peak in peripheral nuclear collisions. The natural 
energy spread of emitted single neutrons[l] being approximately 10% a reso- 
lution of^ < 20% at En= 100 GeVappeared reasonable. 

The limited available space between the RHIC beams at the ZDC location 
imposes the most stringent constraint on the calorimeter design. As can be 
seen from Figure 1, the total width of the calorimeters is only cannot exceed 
10 cm (equal to 1 nuclear interaction length (A/) in tungsten). We designed 
the ZDC's to minimize the loss in energy resolution due to shower leakage, 
which can cause fluctuation in measured shower energy through dependence 
on position of impact and random fluctuations in shower development. 

Finally, the ZDC's are required to withstand a dose of ~ 10^ rad., which is 
the expected exposure during several years of RHIC operation[3]. 
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Fig. 2. Apparent shower profile for different hadron calorimeters. The W and 
Cu parameters are described in the text. A variant with numerical aperture=0.2 is 
shown for comparison. 

3 Simulations 

We simulated shower development, light production and transport in the opti- 
cal components using Geant 3.21 [4] for 2 basic sampling calorimeter designs: 

(1) Pb absorber with scintillator sampling 

(2) Pb, Cu or W absorber, each with undoped fiber optical ribbons in the 
sampling layer 

The ZDC sampling technique which we adopted for this project, is sensitive 
to Cerenkov light produced by charged shower secondaries in a commercial, 
PMMA based [12] communication grade optical fiber with characteristics given 
in a single free parameter was used to match the (wavelength dependent) 
optical fiber attenuation coefficient and photomultiplier quantum efficiency to 
the observed signal from testbeam // 's. Hadronic shower simulation is based 
on Geisha [4] with a low energy cutoff of 0.5 MeV on electrons and photons 
and 1 MeV on hadrons. 

The fiber sampling layers, in all cases, consist of a single ribbon of 0.5mm 
diameter fibers as shown in Fig.5. We chose an orientation of 45° relative 
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Table 1 

Mechanical parameters of the ZDC's 





Absorber 


Space for fibers 


Modules /Layers 


Prototype W-ZDC 


Tungsten 
(100x150x5 mm^) 


1.0 mm 


4(8A/;218Xo) 
27 


Prototype Cu-ZDC 


Copper 

(100x150x10 mm^) 


1.0 mm 


8(7.5A7;79Xo) 
10 


Production ZDC 


Tungsten alloy 
(100x187x5 mm^) 


1.4 mm 


3(5.1A/;149Xo) 
27 



Table 2 

Characteristics of the fiber rib 


3on material. NA = 0.50 




Outer Diam 


Material /Index 


Core 
Cladding 
Surface Prep 


0.45 mm 

0.50 

0.60 


PMMA/1.49 

fluorine doped/ 1.40 

White SiHcone Rubber/ EMA 



to the incident beam direction which roughly coincides with the Cerenkov 
angle of P—1 particles in PMMA. PMMA fibers are readily available with a 
numerical apcrture(NA) of 0.5 (defined as the fractional solid angle which is 
transmitted in the fiber). This aperture corresponds to a maximum angle of 
30". Quartz fibers generally have a smaller aperture. 



In our simulations we studied: 



(1) the effects of transverse shower leakage 

(2) energy resolution dependence on sampling frequency and photostatistics 

(3) dependence on fiber orientation 

Figure 2 illustrates the main advantage of the zdc(Cerenkov) vs. scintillator 
sampling technique. We plot the fraction of the calorimeter signal (photons 
transported in the scintillator or fiber) produced as a function of radius from 
100 GeV protons impacting the calorimeter. A radius of 5 cm (the maximum 
space allowable at the RHIC location) contains 75% of the shower signal in the 
case of a Pb/Scintillator calorimeter with lOmm(Pb) and 2.5mm(Scint) layers. 
The same dimension Pb calorimeter with Cerenkov sampling yields 91% con- 
tainment. In general the Cerenkov technique with Pb absorber achieves a given 
level of containment with a factor of 2 smaller radius than with scintillator. 

Changing from Pb to W absorber yields almost another a factor of 2 reduc- 
tion in containment radius. On the other hand, reducing the fiber numerical 
aperture in the sampling layer produces only a negligible change. 



4 



^ 0.8 

a 

2 0.7 
U 

^ 0.6 

O 0.5 
O 

J3 0.4 
0.3 
0.2 
0.1 






Electrons 




a) 


■ 


P=0.60Mev 
P=0.75Mev 
P=1.00Mev 
P=50.0Mev 




1 \ 














"it. \>_ 



20 



40 



60 



80 



1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 




: 1 

Muons 




^ I 


b) 


— \ ■ 


Beta=0.75 


\ A 


Beta=0.85 




Beta=0.90 


/ X/'^^ * 


Beta=0.99999 






.1/1 1 1 1 1 1 





20 



40 



60 



80 



a 
o 
;-i 
+j 
u 

o 
o 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




L Beta=0.99999 


L c) r 

[ 1 


\ * NA=0.5 
\ A NA=0.4 
\ ▼ NA=0.3 
\ o NA=0.2 

\\ 






^ / ^ 1 




— / •* / / 


\ V * A 


* y Q) 






1 V> Nl 



20 40 60 80 

Angle(Dgr) 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 





20 40 60 80 

Angle(Dgr) 



Fig. 3. Photoelectron yield for 0.5 mm diameter PMMA fibers. 
3.1 Cerenkov light Production and Capture 



Because the optical fibers only transport Cerenkov light emitted nearly aligned 
with the fiber axis, this detector is most sensitive to charged particles which 
cross at approximately 45° to the fiber axis. The lower energy shower compo- 
nent, which is more diffuse, is therefore suppressed. 

This filtering effect is reduced by multiple coulomb scattering of electrons and 
by the increased path length traversed by particles with less than 45° angle 
to the fiber direction. 

In Figure 3 we plot the average photoelectron yields for muons and electrons 
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traversing a ribbon of 0.5mm diameter fibers as a function of incidence angle. 
Comparing Figures 3a and 3b we see that multiple coulomb scattering has 
a significant effect on the response to low energy electrons. The improved 
angular filtering that could be obtained with lower NA fibers (see fig. 3c) is 
largely offset by this effect and the lower light yield which results from small 
aperture. 

Fig. 3d illustrates the effect of "cladding modes" which we suppress by ap- 
plying an ExtraMural Absorber (latex based black paint) over the cladding 
layers of our fibers. 

3.2 Relative response to electrons and hadrons 

Figure 4 shows the calculated response of our zdc modules to electrons and 
protons as a function of beam energy. The response to beam energy (100 GeV) 
protons is a factor of 2 lower than for electrons of the same energy. In this 
sense, our design is extremely non-compensating. This lack of compensation 
is the dominant source of energy resolution of the calorimeter for 100 GeV 
protons since the response changes with fluctuations in the energy fraction 
carried by 7r°'s in the hadronic shower. 

3.3 Linearity 

In our application, where the calorimeter is used to count beam energy neu- 
trons, linearity is not a design consideration. It is clear from figure 4 that , 
while the response to electrons is linear with energy, the hadron response is 
not. Response to neutrons and protons approaches zero at low energies. Also 
shown in figure 4 is the response to muons, which is energy independent up 
to ~ 100 GeV/c at which point radiative energy loss becomes significant. We 
use cosmic ray and beam muons for detector pre-calibration. 

3.4 Energy resolution 

The role of the main components of the energy resolution is illustrated in 
Table 3, where we fit simulated calorimeter response to 50-800 GeV proton 
induced showers to a stochastic plus a constant term. Our results are a poor 
fit to a quadratic sum. 

If we were to increase the sampling frequency from 1/2 to twice our design 
value of 1 per 5mm tungsten absorber there would be a negligible change in 
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Fig. 4. Calculated response of the ZDC to protons and electrons(for both Tungsten 
and Copper modules) and muons (shown for the tungsten modules). 

resolution at 100 GeV. This configuration would reduce the stochastic term 
due to photostatistics from 6 to 3% but leave the dominant resolution term, 
due to non-compensation, unchanged. 



4 Module Construction 



For 10 cm wide modules with 5mm absorber plates, a convenient longitudinal 
segmentation is 1 module per 2 nuclear interaction lengths of absorber. The 
total fiber area matches that of a standard 2" PMT. 
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Table 3 

Calculated energy resolution of ZDC's 



Absorber 


PhEls 


"e/h" 


Stochastic 


Constant 




per lOOGev 


ratio 


term(%) 


term(%) 


W (2.5mm) 


1036 


1.79 


69.6±7.9 


10.1±0.7 


W(5.0mm) 


518 


1.78 


84.6±4.8 


9.1±0.5 


W(lOmm) 


256 


1.78 


92.4±8.2 


8.8±0.6 


Cu(lOmm) 


611 


2.01 


111.7±7.0 


9.3±0.6 


Pb(lOmm) 


422 


1.80 


91.0±8.9 


9.5±0.6 




Fig. 5. Mechanical design of the production Tungsten Modules. Dimensions shown 
are in mm. 

For the prototype W modules we obtained 2.5 mm thick cast plates from a 
Russian manufacturer and bonded them in pairs. For the production mod- 
ules we obtained machined plates of tungsten alloy with threaded mounting 
holes from a US manufacturer [10]. The thickness uniformity of our plates is 
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± 0.1mm. 



The fiber ribbons were wound on a mandrill and then impregnated with a 
low viscosity white silicone rubber glue [9]. The glue covers the active region 
of the fibers (200 mm) and protects the fiber surface in the region of the 
fiber /absorber sandwich. The light guide section of the remaining fibers is 
treated with an extramural absorber to suppress cladding modes in the fiber. 
The fiber ends closest to the PMT are collected into an acrylic compression 
fitting and impregnated with epoxy (Bicron BC-600) . After the epoxy cured 
the fiber bundle was polished using a diamond tipped cutting tool on a milling 
machine. 

The far end of the fibers were rough cut and left untreated. Our optical sim- 
ulations assume no reflection at this end. 

We removed 3 fibers at random from the ribbons in each module and coupled 
them to a single external optical connector for PMT gain monitoring. This 
allowed for stable optical connections of all modules in the calorimeter stack 
to a single hght flasher and therefore rehable tracking of relative PMT gain. 

As illustrated in Fig. 5, the fiber ribbons were trimmed to different lengths 
depending on their positions along the module. Lengths were adjusted to com- 
pensate for the difference in arrival time between the front and back of the 
module. We kept the length of the acrylic fibers to a minimum because we 
were concerned about additional light production in fibers outside the absorber 
region- primarilly due to shower leakage at the top of the calorimeter. 

We selected a 12-stage general purpose PMT (Hamamatsu R329-2)[13] and 
mounted it with a 0.5mm air gap from the fiber bundles. PMT's were se- 
lected for < 6% photocathode non-uniformity over the 39mm diameter area 
corresponding to the fiber bundle size. Linearity of the PMT/ voltage divider 
combination is also an important criterion for this project since the calorime- 
ters will be used to measure up to 40 or so beam energy neutrons in collisions 
of gold ions. 



5 Precalibration 

All modules were tested for relative light yield using cosmic muons incident 
along the beam axis and a standard PMT with cahbrated response. Very httle 
variation (< 10%) was observed among the 24 tungsten modules we eventually 
installed for this project. 

During testbeam operation we also installed trigger counters to select a small 
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contamination of muons which were present in the proton beam. Muons travers- 
ing the fuU calorimeter were used to adjust relative gains of the 4 PMT's in 
the prototype calorimeter. 



6 Testbeam performance 

The calorimeters were mounted on a table with remote x-y (transverse to 
the beam) positioning in the CERN North area, downstream of experiment 
NA49. The beam size and position was defined to ±1 cm using a 1cm square 
scintillation counter directly upstream of the table. 

The main purpose of the beam test was to study the response and resolution 
of the calorimeter as a function of position. 

The 100 and 160 GeV/c protons were selected using a beam Cerenkov counter. 
The beam energy spread was typically 1%. PMT current pulses were integrated 
and digitized using a standard commercial ADC (LRS 2249w). 



7 Results 

Figure 6 shows the measured lineshapes with 100 and 160 GeV incident pro- 
tons. Our energy scale is normalized using the 100 GeV point. Also shown in 
Fig. 6 are the expected distributions based on Geant simulation. The distribu- 
tions are well represented by a gaussian resolution function and the response 
is linear over this limited energy range. 

In Figure 7 we plot the prototype calorimeter response vs. impact position for 
both W and Cu module types. The position scans show essentially uniform 
response to within 1cm of the calorimeter edge along the horizontal direction 
-in good agreement with simulation. In the vertical scans there is an abrupt 
increase in response near the upper edge of the modules. Our simulation re- 
produces this edge effect. It can be traced to shower leakage into the fibers 
above the absorber. 

In order to improve module uniformity in the "beam region" we increased the 
height of the module in our final production design from 10 cm to 13.6 cm. 

The tungsten calorimeter uniformity and energy resolution were essentially 
unchanged when the energy deposited in the 4th module was neglected. Typ- 
ically 1-2% of the energy is seen in this module. The energy resolution at 100 
GeV changes from 17.6 to 19% when it is removed. 
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Fig. 6. Tungsten ZDC response linshapes for 100 and 160 GeV incident protons. 

In Figure 8 we plot the resolution (r.m.s./mean) for both prototypes (using 
the full 8 Aj) at 100 and 160 GeV and compare them to simulation. Part of 
the resolution of this calorimeter is due to the unequal response to electrons 
and photons relative to hadrons. This introduces a limiting resolution due 
to fluctuations in the shower composition (ie vr*^ vs charged vr). Hence the 
non-zero intercept in Fig. 8. 

Our calorimeter is designed to measure beam energy neutrons incident at the 
front face (45° to the fiber direction) . A by-product of the directional response 
of the calorimeter is that it is relatively insensitive to background particles 
from "beam halo" and other sources. In order to demonstrate this suppres- 
sion, we inverted the calorimeter in the testbeam (135° to fiber direction). 
The expected performance is shown in Fig. 9 together with the measured re- 



11 




1 2 3 4 5 

X(cm) 




-4 -2 2 4 



Y(cm) 

Fig. 7. Response map of the calorimeters. 

sponse in the 2 configurations. The small discrepancy in the 135" data could 
be accounted for by a 10% reflection coefficient at the open ends of the fiber. 

1.1 Timing Resolution 

The ZDC signal formation is intrinsically fast. The main spread in transit 
time of light emitted in the fibers results from propagation along the fibers 
and depends on the vertical position of shower particles. A design goal of the 
calorimeters is to achieve better than 200 psec time-of-fiight resolution. Time- 
of-fiight measurement from the interaction to each of the ZDC's is useful 
since it can be used to determine the origin of beam interactions (which is 
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Fig. 8. The ZDC energy resolution versus energy, 
proportional to the time difference). 

In practice, the arrival time of the shower will be measured by digitizing 
the time at which 1 or more PMT signals cross a threshold. We corrected 
the measured time taking into account the signal risetime(2.5 nsec) and the 
measured amplitude. 

Fig. 10a shows a fit to the testbeam data from which we determine the slew 
correction. In fig. 10b we plot the time difference between signals from modules 
1 & 2. In figures 10c, d we plot the time difference of the two module signals 
compared to the beam defining scintillation counters. The resolution of the 
beam defining scintillators is no better than that of the zdc but it is plausible 
, based on figure 10b, that the intrinsic ZDC resolution is Ot < 150ps and it 
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Fig. 9. The ZDC response versus angle. 

is certainly better than 200 psec. With a 150 psec resolution the interaction 
point is measured to an error of 3 cms which is perfectly adequate for our 
application. 

1.2 Radiation tolerance 

PMMA is not a particularly radiation tolerant plastic. It is known to lose 
transparency more readily than Polystyrene, for example. Earlier measure- 
ments on acrylic fibers [11] showed about a factor of 2 decrease in attenuation 
length per 10^ rad of gamma irradiation. Doses at the ZDC location in RHIC 
have been estimated at lOkrad/yr. This is confirmed by dosimetry studies 
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Fig. 10. Time resolution. 



during RHIC commissioning. 



We exposed one of our prototype modules to much higher integrated doses at 
a reactor [14]. Approximately 2/3 of the total dose was due to gamma rays and 
the remainder was due to neutrons The module light output was measured 
using cosmic ray muons before and after 3 successive exposures up to a max- 
imum of 700 krad (7000Gy). The results are plotted in fig. 13, from which we 
conclude that the useful lifetime of the calorimeters will be > 500 krad. 



15 



100 

o 

> 



90 



PSj 80 



70 



1 



1 



60 



T 



200 400 600 800 1000 

DOSE(KRads) 

Fig. 11. Shower response after after radiation exposure at a reactor. 
8 Production Design Choices 



The Copper and Tungsten modules both had adequate performance for our 
application. The tungsten module yields 1-2% better energy resolution at 100 
GeV and slightly better flatness of response over the calorimeter face. We chose 
to proceed with the Tungsten module design primarily because of the 2 x A/ 
modularity and other aspects of the mechanical design. Also the modules are 
more compact, which is an advantage given our limited space. The module 
height was increased in the production design to reduce shower leakage into 
the flber bundles. Also the calorimeter depth was reduced from 8 to 6 x A/ 
since the shorter calorimeter gives essentially identical performance. 
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Fig. 12. Online Zero Degree Calorimeter Energy Distribution obtained during 
RHIC colliding beam operation with beam energies of 65 GeV/nucleon. 

9 Discussion 



Studies of the Cerenkov/fiber sampling technique with electrons [8] and have 
been reported elsewhere. One device [5] was built and operated in a fixed 
target experiment. Here we report first measurements with a 45° design hadron 
calorimeter and the first application of this method at a collider. 

During the first colliding beam operation at RHIC, with Gold Ions accelerated 
to and stored at 65GeV/nucleon beam energy, the ZDC's have been used for 
beam tuning and as a trigger by the RHIC experiments. Figure 12 shows 
an online Energy distribution from one of the calorimeters in the PHENIX 
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experiment. The single neutron peak, which is seen clearly in this distribution, 
has been used to confirm the energy calibration of the calorimeters. 
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